REPORT No. 791 



A THEORETICAL INVESTIGATION OF LONGITUDINAL STABILITY OF AIRPLANES WITH FREE 
CONTROLS INCLUDING EFFECT OF FRICTION IN CONTROL SYSTEM 

By Hakry Gebenbbhg and Leonabd Steenfield 



SUMMARY 

The relation between the devator hinge-momerd parameters 
and the control forces Jor chaiiges injorvxvrd speed and in maneU" 
vers is shovm jor several values of static stability and elevator 
mass balance. 

The stability of the short^eriod oscillations is shovm as a 
series of boundaries giving the limits of the stable region in 
terms of the elevator hinge^moment parameters. The effects 
of static stability J elevator inoment of inertia^ elevator ma^s 
v/nbalance, and airplane density are also considered. Dynamic 
instability is likely to occur if there is mass unbalance of the 
elevator control system combined with a small restoring ten- 
dency {high aerodynamic balance). This instability can be 
prevented by a rearrangement of the unbalancing weights which, 
however y involves an in^ease of the amount of weight necessary. 
It can also be prevented by the addition of viscous friction to 
the elevator control system provided the airplane center of 
gravity is not behind a certain critical position. 

For high values of the density parameter j which correspond 
to high altitudes of flighty the addition of moderate amouTiis of 
viscous friction may be destabUi^ng even when the airplane is 
sialicaUy stable. In this case, increasing the viscous friction 
makes the oscUlalion stable again. The condition in which 
viscous friction causes dynamic instalnlity of a statically stable 
airplane is limited to a definite range of hinge-moment param- 
eters. It is shown that, when viscous friction causes increas- 
ing oscillations f soUd friction wUl produce steady oscillations 
having an amplitude proportional to the amount of friction. 

INTRODUCTION 
Tho effects of aerodynamic balance and mass unbalance 
of the elevator on the djoiamic stabiUty of the airplane are 
discussed in a previous report on control-free stability (refer- 
ence 1) . It was found theoretically ia reference 1 and verified 
in flight (reference 2) that, if the elevator is too closely 
balanced aerodynamically and has a sufficient amount of 
miiss unbalance (which tends to depress the elevator), in- 
creasing oscillations of short period may occur. Other flight 
tests (reference 3) showed, however, that mass unbalance of 
the elevator control system improves the static stability of 
an abplane, that is, increases the slope of the curve of stick 
force against speed in level flight and of the curve of stick 



force against normal acceleration in maneuvei-s. Subse- 
quent work (reference 4) has indicated that a control surface 
with positive floatiag tendency (tendency to float against 
the relative wind), when used as a rudder, is effective in im- 
proving control-free static stabihty. A theoretical analysis 
(reference 5) showed that a rudder having a positive floating 
ratio may, imder the influence of soHd friction in the control 
system, build up steady oscillations of the airplane and 
rudder. These steady oscillations have been observed in 
flight tests (reference 6). These results suggested an investi- 
gation of the behavior of an airplane equipped with an 
elevator having a positive floating tendency. This type of 
elevator was not considered in any of the previous investi- 
gations. 

The purpose of the present report is to make a theoretical 
analysis of the control-free longitudinal stability of an air- 
plane, which takes account of this current trend toward a 
positive floating tendency in control-surface design and 
covers, in general, a much wider range of parameters than 
the investigation of reference 1. These parameters include, 
for the elevator control system, restoring tendency, floatiag 
tendency, mass unbalance (bob weight control), moment of 
inertia, and viscous and soUd friction and, for the airplane, 
density and center-of-gravity position. 

The method of analysis of dynamic stability is based on 
the classical theory of Bryan and Bairstow extended to in- 
clude movements of the controls and their couplings with 
the airplane motions. Friction is treated in the same way 
as in the approximate method of reference 5, in which solid 
friction is replaced by an equivalent viscous friction. 

Before the analysis of dynamic stability is presented, some 
discussion is given of the effect of the various parameters on 
the elevator forces for trim and for acceleration — charac- 
teristics considered important to flying qualtities. The 
stability of the short-period oscillations, with and without 
friction ia the control system, is then considered. The effects 
of weights added to the system to modify the static and dy- 
namic stabflity are discussed. The trends to be expected 
are illustrated by a series of calculations and charts based 
on a typical airplane. The stability of the long-period 
(phugoid) osciUations is not discussed because of its rdative 
unimportance. 
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SYMBOLS 

A„ wing aspect ratio 

At tail aspect ratio 

Ay Bf C, E, F coefficients in stability equation 

b wing span 

Ofi elevator hinge-moment coefficient 

E 



frictional 



(k 

G„ 
c 

Cm 

D 

F 

F, 

Fu 

9 
H 

B, 

H. 

Br 

h=rh,+ht 
4B. 



h = 

I. 



tl=t« — 7%! 

87. 

kr 



hinge-moment 



coefficient 



applied hinge-moment coefficient 
airplane lift coefficient (^^^ 
lift coefficient of tail 

pitching-moment coefficient about airplane 

center of gravity 
wing chord 
elevator chord 

differential operator 

constant term in stabUity equation 
stick force; positive for pull 

(dF \ 
~dn) 

(dF \ 

acceleration of gravity 

hinge moment; positive when tends to de- 
press trailing edge 

mass moment of elevator about its hinge; 
positive when tailheavy 

mass moment of control stick about its pivot; 
positive when stick tends to move forward 

frictional hinge moment 



moment of ioertia of elevator about its 
hinge 

moment of inertia of control stick about its 
pivot 



2i* 



radius of gyration of airplane about I^-axis 



distance between airplane center of gravity 
and elevator hinge 



M 
m 

n 



P 

a 

Stc 

8 

Tin 
t 

AV 

V 
AV 

w 

X 



8e 



T 

e 

X 

/* 

p 



length of control stick 

pitching moment about airplane center]^ of 

gravity 
mass of airplane 

number of cycles required for oscillation to 
damp to half amplitude 

normal acceleration per g of airplane due to 
curvatme of ffight path; accelerometer read- 
ing minus component of gravity force 

period of oscillation, seconds 

dynamic pressure 

elevator area 

tail area 

wing area 

distance in half-chords i2Vt{c) 

time required for oscillation to damp to half 

ampHtude, seconds 
time 



forward velocity 

change ia forward velocity from trimmed value 
weight of airplane 
longitudiual force; positive forward 
distance of center of gravity from aerodyna- 
mic center; positive when center of gravity 
is ahead of aerodynamic center 
normal force; positive downward 
angle of attack 
angle of attack at tail 

deflection of elevator; positive for downward 

motion of trailing edge 
amphtude of elevator oscillation 
angle of downwash 
control gearing i9j5,) 
angle of pitch of airplane 
deflection of control stick; positive for for- 
ward motion of stick 
complex root of stabiUty equation 
real and imaginary parts, respectively, of X 
airplane-density parameter (m/pSJ)) 
mass density of air 
Whenever a, at, B, b, Da, Vd, D5, and D^a are used 

as subscripts, a derivative is indicated. For example, 



Whenever a dot is used above a 



sjonbol, it denotes differentiation v/ith respect to time. 
All angles are measured in radians. 

METHOD OF ANALYSIS 

Four degrees of freedom — forward speed, angle of attack, 
angle of pitch, and elevator deflection — are generally in- 
volved in the problem of control-free stability. To each 
degree of freedom, there- corresponds an equation of equilib- 
rium between inertial and aerodynamic forces or moments. 
By use of wind axes, the four equations become, for level 
flight, 
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XyAV+XaCt 

Zv^V+Zaa 

MyAV+MaOC+Mia+M-a'd 
HyAV+Haa+Hid+H-ad 



+X,d =mAV 

=mV{d—0) 
+Mi$+MiSe+MtS,=mkr^d 

+Hie+HsS.+HiS. =I»(e+5,) +E,{Vd- Ve+L^i) +r[/.(r5.-6») +H.(yd- v'e)] 



(1) 



whidi can be written in nondimensional form as 

(—x^+2A^fiD)u—Xaa =0 

Ci,u+(^^+2A„liD^ a -2A^nDd =0 

In applying equations (1) to dynamic stability, certain approximations may be made. For instance, sbort-period oscil- 
lations (of ttie order of 1 sec) involve negligible changes in forward speed, which, may therefore be neglected in studying the 
shorfc-period oscillations. In fact, the period and damping of these oscillations can be obtained to a high degree of accuracy 
by using only the last three of the equations (1) and setting %=Q. 

Equations (1) then become 

(^-\-2A„yI)^a -2A„iJ)e 



{Oi„^+G„^J)+C„^jy^a +{C„^-2A„yk,W)DB + ((7„,+<7„^/>)5, =0 



(2) 



By setting 



it can be shown (reference 7) that X must be a root of a quartic equation formed by writing 

—2A„iJi 0 



Q 

-^+2A„,x\ 



Gi>t„+ 0„j^^\+ C'«^„X' 



G,„+(G,^^-h)\+G,^J.' G,^+h-(hji,+i{)\ G,,+G,^,\-i,\' 
The resulting stabiUty equation may be written as 



=0 



(3) 



where A, B, O, E, and F axe functions of the stability derivatives. 



The study of the eflFects of different parameters on the 
control-free stability was made by a series of computations 
for an average airplane having the characteristics given 
hereinafter. The current trend toward a positive floating 
tendency in control-surface design suggests the use of Ojt^^ and 
Of,^ as the fundamental variables to be used in expressing 
stability and control characteristics. The results are pre- 
sented as a series of figures that show the relations between 
and Ojif which, with the other derivatives fixed, satisfy 
the conditions for neutral dynamic stabihty and neutral 
static stability. 

A curve for neutral dynamic stability is the boimdary 
dividing the region of increasing oscillations from the region 



of damped oscillations and is obtained from Routh's 
discriminant 

The condition for neutral static stability is that 

The stability equation (3) has four roots. A pair of 
complex roots indicates an oscillatory mode and a real root 
indicates an aperiodic mode. The real part of the complex 
root determines the damping; the imaginary part determines 
the period of the oscillations. More specifically, if there is 
a pair of complex roots 
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the period in seconds is given by 

and the time in seconds to damp to half amplitude is given by 
rp c 0.693 

For an airplane at constant speed, there may be two 
oscillatory modes, there may be only one oscillatory mode, or 
the motion may be entirely aperiodic. In cases in which 
there are oscillatory components, one of the osdllations may 
be poorly damped and even become imstable. 

The average airplane on -which the calculations of this 
report are based is of conventional design. The char- 
acteristics of the airplane are 

- 6 

Jk, L5 

r 1 

sjs^. a 18 

ft 2 

SJSt^.. 0.56 

At 4. 5 

The basic stability derivatives and parameters obtained from 
these airplane characteristics by methods shown in appendix A 
are 

Cl^ 4.3 -16.3 

Cl,^^ — 3.8 a„^, -8.9 

0.486 3.22C*„^ 

-1.64 C„^.— 2a 2 

C^j^ -0-97 e,^^ -10.55a„^ 

Chjii (with no friction) — — 1 

The following parameters of the airplane were varied: 

II airplane-density parameter 

C7m„ control-fixed static-stabihty parameter 

The following parameters of the elevator control system 

were varied: 

Ofl^^ floating tendency 

Cvkj restoring tendency 

Ca^j elevator-damping parameter 

% mtoment-of-inertia parameter of elevator about its 
hinge 

if moment-of-inertia parameter of control stick about 
its hinge 

h mass-moment parameter of elevator control system 

about elevator hinge 
A, mass-moment parameter of elevator alone about its 

hinge 

As has been pointed out, the stability boimdaries were 
plotted, in most cases, in terms of On^^ and O^^ as the varia- 
bles of the coordinate system. In analyzing the effects of 
friction in the control system, Cij and 0*^, were used as the 
plotted variables in some figures whereas C»^^ and Oij were 
used in others. The effect of the other parameters is found 



by varying them one at a time, through a range of values, and 
showing for each parameter a series of stability boimdaries. 
The size of the airplane, wing loading, and altitude are 

combined in the parameter /x, which is '^^j^' A variation 

in PL thus could be due to a variation in size, wing loading, 
or altitude, or any combination of these. The range of values 
of /X covered in the present report and some typical correspond- 
ing values of wing loading, altitude, and size are given in the 
following table: 





Wing 
loading 
Clb/«q ft) 


Altitude (ft) 


Mean 
wing 
ohoru 
(ft) 


4.17 
12.6 
37.6 


40 
40 
40 




21 
7 
7 




33,000-- 





The range of and the corresponding center-of-gravity 
positions are as follows: 





(fraction of 




mean wing 




chord) 


-0.232 


0.05 


0 


0 


.233 


-.06 



The ranges of values of the other parameters, for a small 
airplane (chord, 7 ft), are as follows: 





Moment of 




inertia of ele- 


f« and U 


vator control 




system 
(sfug-fti) 


0 


0 


2 


1.8 


4 


3.2 



h 


stick force (lb) 


At sea level 


At 33,000 ft 


0 
10 


0 
37 


0 

12H 





dH/dS 




V 




Gb-ft/lOO mph/deg/sco) 




At sea level 


At 33,000 a 


—10 


6.86 


1.06 


—100 


5S.5 


19.6 



STATIC STABILITY AND RELATION TO CONTROL FORCES 

The connection between the static stability and the 
airplane and control parameters is established to assist 
in the interpretation of the results obtained hereinafter. 
Equations (1) can be applied to static stability by setting 
all terms containing D and equal to zero and solving the 
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resulting equations simultaneously - for the variation in 
forward speed with, an apphed elevator hinge moment. 
For level flight, 6 is also zero and the resulting equations are 

Or.U+^ a =0 

Solving gives 

The variation of stick force with fractional change iu forward 
speed is 

If effects of sUpstream on the tail are neglected, (7«y=0. 



As shown in appendix A, ^a«=~2Xm' Inserting these 
values in the expression for Oj^Ju shows that Fu is iadependent 
of forward speed. 

The variation of control force with normal acceleration in 
a steady pull-up, with no change assumed in forward speed 
(see reference 8), can be foxmd from equations (2) by equating 
to zero aU terms containing D except Dd. This procedure 
imphes that the normal acceleration is due entirely to 
cmrvature of the flight path Dd, The equations become, for 
an applied hinge moment, 

(J 

oc-2A„fjLDd =0 

from which 



D9- 



If the normal acceleration is Tig, 



and 

DO 

These formulas for Fu and Fn are equivalent to equation (1) 
of reference 9 and equations (27) and (28) of reference 10. 

The formulas indicate that the stick-force gradients F« 
and Fn are dependent on most of the aforementioned airplane 
and elevator parameters. Figures 1 to 5 show the variation 
of these stick-force gradients with the parameters Cii^^, On^, 
Omai ^* The gradients are independent of speed, 

although only within the limits of the assumptions made in 
the analysis, namely, neglect of power and of compressibility 
effects. The gradient Fu can be used to get the stick force 
for only a small change in forward speed because the stick 
force is not directly proportional to the change in speed. 
The stick force in a steady pull-up Fn, however, is propor- 
tional to the normal acceleration provided the control de- 
flection is not so great that the basic assumption of linearity 
is violated. 

The line Fu—0 is the boundary for true static stability — 
that is, jP'u=0 is the condition for zero variation in stick 
force with forward speed in steady flight. This condition 
is the same as that obtained by setting F=Q, where F is 
the constant term of the sixtii-order stability equation 
obtained fi*om equations (1). On subsequent figiures it is 



called the divergence boundary. The line Fn=0 is the 
boundary for apparent static (or maneuvering) stability and 
is the condition for zero variation in stick force in a steady 
pull-up. This condition for Fn=0 is obtained by setting 
^^=0 in the approxioiate stabihty equation (equations (3)), 
which is for three degrees of freedom (a, Dd, and 5,). On 
the unstable side of ^"^=0, a slow divergence occiub that is 
noticed by the pilot as an unstable variation of stick force 
with forward speed. The stick force due to normal accelera- 
tion in a pxdl-up is stable, however, unless the conditions are 
such that the airplane is operating on the imstable side of 
Fn=0. 

Figures 1 to 5 indicate that the parameters have the same 
effect on jF„ and Fn except that the altitude affects only i^„. 
They show that the stick-force gradients on an airplane of 
given tail size and center-of-gravity position may be in- 
creased by making the floating tendency C*^^ more positive 

or by mass unbalancing the elevator control system to 
depress the elevator (make it tailheavy). The effect of the 
restoring tendency C^j on the stick-force gradients depends 
on the relative position of the center of gravity and the 
aerodynamic center. If the center of gravity is ahead of 
the aerodynamic center (airplane stable with controls fixed), 
increasing the magnitude of Oj,^ increases the stick-force 
gradients. If the center of gravity is behind the aerody- 
namic center, this effect on Fu is reversed; the effect on Fn 
is not reversed, however, until the center of gravity is well 
behind the aerodynamic center (in this case, about 0.05c at 
sea level and 0.02c at 30,000 feet). If (7*^=0, the stick 
forces are independent of the position of the airplane center 
of gravity. 
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Increase in altitude wdll either increase or decrease 
depending on tlie Mnge-moment parameters. The solid 
line in figure 5 is the locus of values of C^^^ and Ch^ for which 

Fn is independent of altitude. For points to the left of this 
line, Fn decreases with altitude; for poiats to the right of 
this line, Fn increases with altitude. This line is determined 
by the relation 

which, for the case of figure 5, becomes 

Another method of increasing the stick-force gradient in 
level flight Fu consists in applying a constant hinge moment 
to the elevator by means of a spring or bungee. The effect 
of the spring on the gradient Fu is due to the derivative 
Ghu which depends ia the same way on the constant hiuge 
moment, whether it is caused by a weight or by a spring. A 
bungee, which tends to depress the elevator, will therefore 
increase the stick-force gradient in level flight Fu- The 
effect of the bungee on the stick-force gradient in accelerated 
flight Fn will be zero because its action depends solely on 
changes in forward speed. Its effect on the short-period 
oscillations will be zero for the same reason. 

DYNAMIC STABILITY 
NO FBICnON IN CONTROL SYSTEM 

The stabihty of the short-period osciQations without 
friction is shown in figures 6 to 11, which also show the 
boundaries for true static stabihty (divergence boundaries). 
Figure 6 is an example of a more nearly complete presentation 
of the st-ability data than subsequent figures because it 
shows the variation of damping and period of oscillation 
with the hinge-moment paiameters Gn^^ and 0*, for certain 

fixed values of the other parameters. The dampiag, which 
is proportional to $, increases with the magnitude of Oh^ 

The period, proportional to decreases as Ci^^ uicreases. 

Another way of presenting this additional stability data is 
shown ia figure 7, which gives the number of cycles the oscil- 
lation performs before it damps to half ampHtude. It is 
clear from figure 7 that the oscUlation is very weU damped 
unless the restoring tendency is close to zero. In this parti- 
cular case, only one oscillatory mode exists. Inasmuch as 
there are only three roots ia this case (because ts and ii=0), 
the other root is always real and is of no particular signifi- 
cance in dynamic stabUity. In cases ia which an additional 
oscillatory mode exists, it is always stable. 

The effect of center-of-gravity position on the stability of 
the short-period oscillations is shown ia figure 8. The shift 
in the dynamic-stabihty boimdaxy, for the comparatively 
large shift ia center of gravity shown, is practically negligible. 



Many of the subsequent figures, ia which zero static stability 
is assumed to facilitate computation, therefore are vaUd for 
airplanes having a margin of static stabihty. 

The effect of moment of iaertia of the elevator control 
system on the dynamic stabihty is shown in figure 9, which 
gives typical values of the moment of iaertia. The effect is 
shghtly destabilizing especially for high values of C*^^. The 

destabUiziDg effect of the moment of inertia of tlie elevator 
is greater than that of the moment of inertia of the control 
stick. Because the accuracy gaiaed by iacludiog moment 
of inertia is small compared with the saving in labor gained 
by neglectiag it, moment of inertia of the elevator control 
system was set equal to zero in the subsequent calculation. 
As a result, the stabihty equation becomes a cubic and sub- 
sequent figures are somewhat imconservative. 

The effect of density parameter /* on the dynamic stabihty 
is shown in figure 10. Increase of m has a shght destabilizing 
effect. 

As has been shown, mass imbalanco of the elevator control 
system improves the static stabUity (fig. 4). The effect on 
dynamic stabihty is imfavorable, however, as sho^vn in 
figiu'e 11. The value of mass unbalance shown corresponds 
to a bobweight that is located at the airplane center of 
gravity and requires a balanciag puU of 37 pounds on the 
control stick of a pursuit airplane at sea level. Increasiag 
oscillations occur if the aerodynamic balance is too high 
(low magnitudes of Ch^, especially for negative values of (7^^. 

The effect of the location of the bobweight is shown in 
%ure 12. Each curve represents a different arrangement 
of bobweights and aU arrangements give the same stick 
force. The sohd line corresponds to a weight at the air- 
plane center of gravity (for practical purposes, at the control 
stick). The short-dash line corresponds to a weight at the 
elevator. The long-dash line corresponds to two weights — 
one at the elevator, which tends to make it noseheavy; the 
other at the control stick, which gives the elevator a suffi- 
ciently powerful tailheavy moment that the resultant stick 
force is the same as with the single weight. In the particu- 
lar case represented, the noseheavy moment due to the 
weight at the elevator is equal to the tailheavy moment due 
to both weights. Moving the single weight fixDm the con- 
trol stick to the elevator has a large destabilizing effect. 
Overbalanciag the elevator while the stick force is kept 
constant has a considerable stabilizing effect. This method 
of preventiag iostabUity has the disadvantage, however, of 
increasing the total amoimt of unbalanciag weight required. 
In the case shown, the weight is iacreased to three times its 
original size. Another disadvantage is the rearward shift in 
center of gravity of the whole airplane due to additional 
weight at the elevator. (See airplane parametei-s given in 
"Method of Analysis.") The destabilizing effect of the ia- 
creased moment of iaertia associated with the added weights 
was found to be very small, especially for negative floating 
tendency. 
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EFFECT OF VISCOUS FHICTION IN CONTKOL SYSTEM 

In the pj'eceding section, a constant value of the elevator- 
damping parameter Ghj^^ was assumed. This value was due 
only to aerodynamic damping. The effects of viscous fric- 
tion in the elevator control system are obtained by consider- 
ing Oa^j as an additional variable. This variable can be 
introduced, as in the preceding section, by showing a series 
of boundaries in the Oh^Oj,^ plane for various values of 

Ofij^y The general nature of the effect of friction is shown 
first, however, by presenting boundaries in the Gbfihj)^ 
plane with Oj,^^ constant and some other parameter varied. 

This method of presenting stability boundaries makes it 
easier to show the effects of other parameters such as air- 
plane center-of-gravity position and density when friction 
is introduced. 

The effect of viscous .friction on the dynamic stability, 
for various conditions, is shown in figures 13 and 14 for 
^=12.5 and j[t=37.5, respectively. Figures 13(a) and 14(a) 
refer to the mass-balanced elevator control system; figures 
13(b) and 14(b) refer to the tailheavy elevator control 
system considered in the preceding section. It is shown that, 
if the airplane center of gravity is ahead of a certain point,^ 
the instability caused by the unbalanced elevator can be 
removed by adding viscous friction to the control system. 
This critical center-of-gravity position is behind the aero- 
dynamic center, and its distance from the aerodynamic 
center decreases as the density parameter /x increases. 
When the center of gravity is behind this critical position, 
viscous friction has a destabilising effect. These opposite 
effects of viscous friction are shown in the Of,^^ O?,^ plane in 
figures 16 and 16. "When the center of gravity is slightly 
ahead of this critical position, the effect of viscous friction 
depends on its magnitude and also on the value of G),y The 
addition of a small amount of viscous friction is destabi- 
lizing but larger amoimts are stabilizing. If the aerodynamic 
balance is sufficiently high (Ot^^O) and the viscous friction 
lies in a certain range, increasing oscillations will occur. 
In figure 14(b), for example, if aja.o.=— 0.01c and Cii3=— 0.05, 
the oscillations will be unstable when the elevator-damping 
parameter is in the range from —2.5 to —76. If is 
more negative than —0.086, no amoimt of elevator damp- 
ing can cause increasing oscillations. As the center of 
gravity moves forward, the destabilizing effect of elevator 
damping becomes less and fiLnally disappears. 

The effect of the density parameter /i can be seen by com- 
paring figures 13 and 14. The critical center-of-gravity 
position approaches the aerodynamic center as /i increases. 

di 

» Stace this report was written, this point has been found to be where sometimes called 
the stick-Oied maneuver point. 



When /x=12.5, elevator damping always has a stabilizing 
effect provided Xa,c, is positive. When ju=37.5, elevator 
damping may be destabilizing over a small range of Ca^j 
and Cb^ even when Xa,c, is positive (0.05c). 

When the center of gravity is slightly ahead of the afore- 
mentioned critical position (which is behind the aerodynamic 
center), the conditions imder which elevator damping may 
cause dynamic instability may be advantageously repre- 
sented in the Ojt^Oj,^ plane. If a series of stabiUty boundaries 

are drawn in that plane for various values of elevator damp- 
ing, they will all be confined to a region boimded by a line 
that will be called the boundary of complete damping. An 
illustration of two methods of constructing this boundary is 
given in figure 17. If a series of boimdaries in the Oh^G^ 

plane are drawn for various values of the damping, the com- 
mon tangent of all these curves is the boundary for complete 
damping. This boundary can also be drawn by plotting the 
minimum values of (7^, obtained from plots of the type shown 
in figures 13 and 14 against corresponding values of Oji^^, 
The region in the Ox^Oj,^ plane between the boundaries for 

complete damping and increasing oscillations is the region 
where the addition of viscous friction to the elevator control 
system may cause dynamic instability. 

That a boimdary for complete damping cannot bo shown 
for /i=12.5 if the airplane is statically neutral or stable (xa.c. 
is zero or positive) may be seen from figure 13. It is possible 
however, to show a boundary for complete damping imder 
these conditions for ^=37. 5. Figrn-e 18 shows these boimd- 
aries for Xa,c=^ aJid for the critical value Xa.c.=— 0.017c, 
for both a mass-balanced elevator and a mass-imbalanced 
elevator. The boundaries for increasing osciQations and 
divergence are also shown. For the case of the mass-balanced 
elevator (A=0), the boundary for complete damping is a 
straight line through the origin and therefore corresponds 
to a fixed ratio of the floating and restoring tendencies, 
or floating ratio. Elevator mass imbalance decreases the 
region of complete damping. 

EFFECT OF SOUD FRICTION IN ELEVATOR CONTROL SYSTEM 

The boundary for complete damping has an important 
bearing on the effect of solid friction on dynamic stabihty. 
In order to calculate this effect, the soHd friction is replaced 
by an equivalent viscous friction that would dissipate energy 
at the same rate. This condition gives an equivalent 



for a sinusoidal motion of the elevator with amplitude 5 and 
angular frequency 17. 
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Restoring tendency, 
FiQUBE 17w— Method of constxactlnfi boundary for complete damping. 



It can be shown that if viscous friction is destabilizing, as 
in figures 16 to 18, solid friction may lead to steady osciQa- 
tions having an amplitude proportional to the amount of 
friction. Suppose an oscillation is started with amplitude 
and frequency which result in an equivalent (7*^ that would 
cause increasing oscillations. Let this condition be repre- 
sented by point 2 in figure 19. The amplitude of the oscil- 
lations would then increase until the equivalent Ca^^ de- 
creased to the value that would result in neutrally damped 
oscillations, represented by point 3 in figure 19. If the 
ioitial ampKtude is so low that the equivalent viscous fric- 
tion is in the stable region, as shown by poiat 1, the oscilla- 
tions will die out completely. If the initial amphtude is so 
high that the oscillations are stable, represented by point 4, 
the amplitude will decrease until it reaches a constant value, 
when the equivalent Ca^ is again represented by point 3. 
The value of 0*^^ at point 3 then determines the amphtude 
of the steady oscfllations. By solving formula (4) for the 
amplitude of the steady oscillation is obtained as 



where i; and are the values at point 3. 



This formula 



shows that the amphtude is proportional to the amount of 
soKd friction. 

The foregoing analysis shows that the region in the Oh^Gh^ 

plane between the boundary for inci*easing oscillations and 
the boimdary for complete damping is the region where 
steady oscillations may occur because of the influence of 
soKd friction in the control system. All the remarks in the 
preceding section concerning the boundary for complete 
damping with viscous friction consequently apply to the 
boundary for steady oscillations with solid friction, inasmuch 
as these two boundaries are the same, within the limits of 
the assmnptions involved in the use of the concept of equi- 
valent viscous frictioD. Steady oscillations due to solid 
friction will not occur on a statically neutral or stable air- 
plane, for instance, unless n is very large (corresponds to a 
high altitude). Even in that case, the possibility of steady 
osciQations exists only for a comparatively small range of 
floating ratios. If the airplane is statically unstable by a 
sufficient amount, however, steady oscillation may exist 
over the entire range of floating ratio. 

Some calculations of the amphtude of the steady oscilla- 
tions, expressed in terms of normal acceleration per unit 
frictional force as felt at the control stick, were made by the 
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Restoring tendency, 

(a) A«10. 

(b) A-0. 

FiouBZ IS^Botmdaries for complete damping, inoreasiiig oscfllatUms, and divergence. 



metliod of appendix B . The results are presented in figure 20, 
which shows lines of constant amphtude in the Ob^Cn^ 

plane for an airplane with the center of gravity at the critical 
position referred to in the preceding section. Steady oscilla- * 
tions win therefor occur throughout the entire region where 
stability exists in the absence of friction. The amphtude 
is smallest for high values of 0^^^ combined with high 

values of Ohy 

CONCLUDING REMARKS 

The stick-free static stabihty of a conventional ahplane 
may be improved by making the elevator floating tendency 
more positive or by mass-unbalancing the elevator control 
system to make the elevator tailheavy. Increasing the 
restoring tendency also has a favorable effect provided the 
airplane is stable with stick fixed. If the restoring tendency 
is zero, the stick-free static stability is independent of air- 
plane center-of-gravity position. 

The dynamic stabihty with fiictionless controls depends 
chiefly on the restoring tendency O^^ and on the mass balance 
of the elevator control system. If the elevator control 
system is mass unbalanced (elevator made tailheavy) by an 
offset weight at the control stick and if the restoring tendency 
is too low, increasing short-period oscillations may result. 
This condition can be remedied by the use of two additional 
weights — one at the elevator making it noseheavy, the other 
at the control stick making the elevator suflSciently taUheavy 



that the combined effect gives the elevator the desired amount 
of taflheaviness. 

The addition of viscous friction to the control system will 
prevent dynEmoic instability provided the airplane center of 
gravity is forward of a critical position which is behind the 
aerodynamic center and approaches it as the value of the 
density parameter fx increases. If the airplane center of 
gravity is behind this critical position, viscous friction wiU 
have a destabilizing effect, no matter what the hinge- 
moment parameters are. If the center of gravity is shghtly 
ahead of the critical position, viscous friction may be de- 
stabilizing for a limited range of values of viscous friction and 
the hinge-moment parameters, A low restoring tendency 
and a high positive floating tendency will tend to cause this 
dynamic iostability. When fi is very lai^e (high altitude), 
this condition of steady oscillations can occur even if the 
center of gravity is ahead of the aerodynamic center. 

The presence of sohd friction in the control system may 
cause short-peiiod steady oscillations under the conditions 
for which viscous friction is destabilizing. The amphtude of 
the oscillations is proportional to the amount of friction 
present. 
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FiauBE 20.— Variation of maiimnm nonnal acceloratioa In the steady oscillation with elevator 
hinge-moment parameters, a:..*.--0^ ;*-37^; A«-IO; f,-2; £,«0. 



Derivative (\. — ^The total hinge moment acting on the 
elevator may be expressed as 

At trim H—Q, therefore, 



APPENDIX A 

EVALUATION OF STABILITr DERIVATIVES 
where 



Hjy hpS^^ 
£±0= — — ' 



n ^ 



Using I pV^0iJ3a=mg gives 



hcgpOLSv, _ Wl 
2mg 2Awfi 



Derivative Om^* — ^The parameter (7«^ may be obtained 
from wind-timnel measuroments or, if the position of the 
aerodynamic center of the complete airplane is known, may 
be calculated by the formula 

Derivatives O^^^ and — The derivatives (7„^^ and 
Omj^'i^ arise because of the lag between the change in angle 
of attack at the wing and the corresponding downwash at 
the tail. It is assumed that the downwash at any instant t 

depends on the angle of attack at the instant t^^f the differ- 
once being the time required for the air to move from the 
wing to the tail. If q:=/(Oi relation may be expressed 
as 



Now, 



Hence, 



/a- At) - Atf (f) (0 - 



c=€a a^Ata+^^^2^a— • • • 
or, because Da and a="^ 

and, because a,=ck:— €, 
and 

I * 

The part of the pitching-moment coefficient contributed 
by the taal is 

The lag ejffectively introduces derivatives C«^^, Omj^^y . • . - 
The first two of these derivatives are 



and 



n — — ^ ^ 
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Derivatives C^^, ^^*2ja> ^^i^a — '^^^ derivatdvea Oia, 
Cft^^, and Ghi^^ may be obtained from 



which gives 



Derivative — ^The pitching moment due to pitching 
is made up of parts] due to propeller, wing, fuselage, and 
horizontal tail. The contribution of the tail is by far the 
lai^est and can easily be calculated. 

If the airplane is rotating with angular velocity dj the 

Q 

increase in angle of attack at the tail is yy which results 
in an increased lift on the tail given (in coefl&cient form) by 

The resultant pitching-moment coefficient is 



ri rt -"ft 'tl n T " ' 



C S„ 



V C 8^ 



and expressing Q bs> — Dd and as h gives 



The contribution of the wing depends on the assumed 
axis of rotation (center of gravity) but a fair average value 
win be obtained by assuming that the center of gravity is 
at the wing quarter-chord point. This assumption gives a 
value 

The total pitching-moment coefficient due to pitching 
therefore is 



11^ 



Derivative Gm^* — ^The derivative Gm^ may be measured 
directly in a wind tunnel or may be computed from \vind- 
tunnel data on the value of G^^^ for the horizontal tail by 

means of the formula 

S~~^'^'«2S« 

Derivative Cm^. — ^The derivative C^^^ may be computed 
from 

where (^^^ (^^) obtained from figure 1 

of reference 11, which is based on thin-wing potential-flow 
theoiy. 

Derivative (7a^. — ^The derivative G^j^g is given by 

In the absence of viscous friction in the elevator control 
system, the value of (7*^, may be computed from 

where (^^^ and (^^^) obtained from figure 1 

of reference 11. 

If a dashpot, which has a damping constant of K pounds 
per foot per second and moves a distance of Q feet per 
radian of elevator deflection, is inserted in the control 
system, 



(A2) 



The total value of Gh^ is the sum of equations (Al) and 
(A2). 

Derivatives C^^^ and G^y — The derivatives Gj,^^ and Gh^ 
can be calculated by thin-wing-section theory but the 
results are of doubtful accuracy because of three-dimensional 
and boundary-layer effects. It is therefore best to obtain 
these derivatives from wind-timnel tests. 
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CALCULATION OP NORMAL ACCELERATION DUE TO OSCILLATING ELEVATOR 
Tho normal acceleration of the airplane, which is equal to D{a — 6) in nondimensional units, can be calculated from 



2 I _£!!|_5«_24„^(7„^_,«(44„VV-^7«A2^«M)+t'; {^-lA^iiC„^,^-2A„vkf-f-2A^i.O. 



(Bl) 



where -q is the angular frequency of the elevator. 



The fraction in equation (Bl) can be reduced to an ordinary complex number and the modulus of this nimiber is the 
maximum amplitude of the steady oscillation. The value can be converted to physical units by the formxila 

Normal acceleration per g _ 41^ D(a—d) 4 
Stick friction in pounds ^ pS^^g 5 ^^^Adj 

where (7*^^ is the value of elevator damping required for the condition of neutral dynamic stability. 
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